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Increasing Accuracy in Computed Inviscid Boundary Conditions 

Accuracy is increased through use of higher-order time derivatives. 

John H. Glenn Research Center, Cleveland, Ohio 


A technique has been devised to in- 
crease the accuracy of computational sim- 
ulations of flows of inviscid fluids by in- 
creasing the accuracy with which surface 
boundary conditions are represented. This 
technique is expected to be especially ben- 
eficial for computational aeroacoustics, 
wherein it enables proper accounting, not 
only for acoustic waves, but also for vortic- 
ity and entropy waves, at surfaces. 

Heretofore, inviscid nonlinear surface 
boundary conditions have been limited to 
third-order accuracy in time for stationary 
surfaces and to first-order accuracy in 
time for moving surfaces. For steady-state 
calculations, it may be possible to achieve 
higher accuracy in space, but high accu- 
racy in time is needed for efficient simu- 
lation of multiscale unsteady flow phe- 
nomena. The present technique is the 
first surface treatment that provides the 
needed high accuracy through proper ac- 
counting of higher-order time derivatives. 

The present technique is founded on a 
method known in art as the Hermitian 
modified solution approximation (MESA) 


scheme. This is because high time accu- 
racy at a surface depends upon, among 
other things, correction of the spatial 
cross-derivatives of flow variables, and 
many of these cross-derivatives are in- 
cluded explicitly on the computational 
grid in the MESA scheme. (Alternatively, 
a related method other than the MESA 
scheme could be used, as long as the 
method involves consistent application 
of the effects of the cross-derivatives.) 

While the mathematical derivation of 
the present technique is too lengthy 
and complex to fit within the space 
available for this article, the technique 
itself can be characterized in relatively 
simple terms: The technique involves 
correction of surface-normal spatial 
pressure derivatives at a boundary sur- 
face to satisfy the governing equations 
and the boundary conditions and 
thereby achieve arbitrarily high orders 
of time accuracy in special cases. The 
boundary conditions can now include a 
potentially infinite number of time de- 
rivatives of surface-normal velocity (con- 


sistent with no flow through the bound- 
ary) up to arbitrarily high order. The 
corrections for the first-order spatial de- 
rivatives of pressure are calculated by 
use of the first-order time derivative ve- 
locity. The corrected first-order spatial 
derivatives are used to calculate the sec- 
ond-order time derivatives of velocity, 
which, in turn, are used to calculate the 
corrections for the second-order pres- 
sure derivatives. The process as de- 
scribed is repeated, progressing 
through increasing orders of deriva- 
tives, until the desired accuracy is at- 
tained. 

This work was done by Roger Dyson of Glenn 
Research Center and Ray Hixon of the Insti- 
tute for Computational Mechanics in Propul- 
sion. Further information is contained in a TSP 
(seepage 1). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Commercial 
Technology Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21000 Brookpark Road, Cleveland, 
Ohio 44135. Refer to LEW-17500-1. 


Higher-Order Finite Elements for Computing Thermal Radiation 

Computationally efficient methods yield close approximations of exact solutions. 

Langley Research Center, Hampton, Virginia 


Two variants of the finite-element 
method have been developed for use in 
computational simulations of radiative 
transfers of heat among diffuse gray sur- 
faces. Both variants involve the use of 
higher-order finite elements, across 
which temperatures and radiative quan- 
tities are assumed to vary according to 
certain approximations. In this and 
other applications, higher-order finite 
elements are used to increase (relative 
to classical finite elements, which are as- 
sumed to be isothermal) the accuracies 
of final numerical results without having 
to refine computational meshes exces- 
sively and thereby incur excessive com- 
putation times. 

One of the variants is termed the radi- 
ation sub-element (RSE) method, which, 


itself, is subject to a number of variations. 
This is the simplest and most straightfor- 
ward approach to representation of spa- 
tially variable surface radiation. Any com- 
puter code that, heretofore, could model 
surface-to-surface radiation can incorpo- 
rate the RSE method without major 
modifications. 

tn the basic form of the RSE method, 
each finite element selected for use in 
computing radiative heat transfer is con- 
sidered to be a parent element and is di- 
vided into sub-elements for the purpose 
of solving the surface-to-surface radiation- 
exchange problem. The sub-elements are 
then treated as classical finite elements; 
that is, they are assumed to be isothermal, 
and their view factors and absorbed heat 
fluxes are calculated accordingly. The 


heat fluxes absorbed by the sub-elements 
are then transferred back to the parent el- 
ement to obtain a radiative heat flux that 
varies spatially across the parent element. 
Variants of the RSE method involve the 
use of polynomials to interpolate and/ or 
extrapolate to approximate spatial varia- 
tions of physical quantities. 

The other variant of the finite-element 
method is termed the integration 
method (fM). Unlike in the RSE meth- 
ods, the parent finite elements are not 
subdivided into smaller elements, and 
neither isothermality nor other unrealis- 
tic physical conditions are assumed. In- 
stead, the equations of radiative heat 
transfer are integrated numerically over 
the parent finite elements by use of a 
computationally efficient Gaussian inte- 
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Temperature Versus Position along a plate was computed in a test case exactly and by two variants of 
the finite-element method. 


gration scheme. In this scheme, the radi- 
ant heat transfer is computed at discrete 
points on each surface in the radiation 
exchange. These points corresponding to 
the Gauss points are used in evaluating 
the element matrices. 

The IM is implemented in the follow- 
ing iterative procedure: 

1. Initialize unknowns (temperatures and 
radiative heat fluxes) . 

2. Calculate differential form factors be- 
tween Gauss points on all elements. 


3. Calculate radiative heat fluxes at 
Gauss points on all elements. 

4. Integrate radiative heat fluxes to ob- 
tain a radiative-heating load vector. 

5. Solve for and update temperatures. 

6. Examine the results for convergence. 
If results have not converged to within 
acceptably narrow margins, return to 
step 3. 

In a test problem, an upper plate of 
length L was assumed to be maintained 
isothermal at a temperature of 1,000 R 


(-556 K) and a lower plate of the same 
length was assumed to be placed at the 
same horizontal position at a distance 0.1L 
below the upper plate and allowed to come 
to thermal equilibrium. The figure depicts 
the temperature of the lower plate versus 
position along the plate as calculated by 
the RSE method with cubic-spline interpo- 
lation and by the IM with 16 Gauss points. 
Also shown is the exact solution. Both 
methods show reasonably good agreement 
with the exact solution, with the integra- 
tion method nearly indistinguishable from 
the exact solution over most of the plate. 
In general, the integration method proved 
to be more accurate with respect to spatial 
variation; however, it was also more cosdy 
(longer run times). Both methods cap- 
tured temporal variations equally well. 
These results indicate that the RSE 
method is preferred for efficient analyses 
in which temperature variations are mainly 
temporal, while the integration method is 
reserved for analyses requiring very accu- 
rate resolution of spatial gradients. 

This work was done by Dana C. Gould 
of Langley Research Center. For fur- 
ther information, access the Technical 
Support Package (TSP) free on-line at 
www.techbriefs.com/tsp under the Informa- 
tion Sciences category. 
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